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ABSTRACT 

We show that the mean phase of waves propagating aU the way from the far side of the Sun to the 
front side, as measured by seismic holography, varies with time. The change is highly anticorrelated 
with solar cycle activity and is consistent with other recent results on the variation of the seismic 
radius of the Sun. The phase change that we observe corresponds to a few kilometers difference in 
the seismic solar radius from solar maximum to solar minimum in agreement with inferrences from 
global helioseismology studies. 

Subject headings: Sun: solar radius, Sun: local helioseismology. Sun: seismic holography 



1. INTRODUCTION 

Accurate knowledge of the solar radius is a key fac- 
tor in structure inversions from solar-oscillation frequen- 
cies. Although several values of the photometric radius 
have been reported by different authors, the accepted 
radius for most inversions has been the revised one by 
iBrown fc Christensen-Dalsgaard I (|1998f ). Measurements 
of temporal variation of the photometric solar radius 
are largely controversial. Different techniques have been 
used to determine the variations and the results differ, 
ranging from no variation, to rapid variations, to solar cy- 
cle variations. iKuhn et al. I ()2004[ ) ascribed the apparent 
disparities in the measured optical solar diameter from 
ground-based instruments to the changing terrestrial at- 
mosphere. They suggest that a solar cycle variation of 
the Earth's atmosphere is responsible for the correla- 
tion/anticorrelation between the observable radius and 
t he solar activi t y cycle . 

ISchou et al. I (jl997f ) introduced the concept of the 
seismic radius of the Sun. They obtained the so- 
lar radius from /-mode frequencies and found a dis- 
crepancy with the photospheric standard value of 
-300 km. iTakata fc Cough I (|2001h confirmed the value 
of the seismic rad i us fro m p-mode analysis. Recently, 
iHaberreiter et al~l (120081 ) explained the discrepancies be- 
tween the photospheric optical radius and the seismic ra- 
dius as arising from the difference between the height at 
disk center and the inflection point of the intensity profile 
at the limb. 

Since the seismic radius is inferred from helioseismol- 
ogy inversions, it does not necessarily agree with the 
physical radius. Changes in the stratification of the Sun, 
such as variations of the sound speed or in the subsur- 
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face superadiabatic layers, influence its determination. 
Hence, variations of the seismic radius are better inter- 
preted as a measurement of changes in the cavity in which 
the waves propagate, rather than a physical change of 
the solar radius. Moreover, localized perturbations due 
to magnetic field concentration can also play a role in 

the determination of the seismic rad ius. 

Recent results from,Lefebvre fc KosovichevI (|2005[ ) us- 
ing /-modes and lKholikov fc Hill I ( 2008f ) using very \ow-l 
p-modes, found a change of the seismic solar radius of a 
few kilometers that is anticorrelated with the solar mag- 
netic cycle. Variations of the solar seismic radius could 
have large implications, from explaining the solar cycle 
variations of the observed oscillation frequencies to af- 
fecting the seismically-inferred surface maps of far-side 
activity. 

Using seismic holography, a local-helioseismology tech- 
nique that is here applied to waves that travel from the 
far side to the front side, we find a variation of the seis- 
mic radius with the solar cycle that qualitatively agrees 
with the global results. The use of local-helioseismology 
techniques for this type of study allows us to discrimi- 
nate between the variation being a global phenomenon or 
a localized one associated with surface magnetic activity. 

2. DATA ANALYSIS 



Seis mic imaging was proposed by iLindsev fc Brauni 
(|1990f ) as a technique to map active regions on the far 
hemisphere of the Sun. Since then, the method has 
also been used suc cessfully to explore the subphotospher e 
of active regions (jChou ifeoOCt iBraun fc Lindsev II2000I) 
investigate seismic sources from flares ()Donea et al. I 
l2006t) and characterize flows in the solar subphotosphere 
(|Braun. Birch fc Lindsev "2004^. Phase-correlation seis- 
mic holography applied to high quality data from the 
Michelson Doppler Imager (MDI) aboard the Solar He- 
liospheric Observatory rendered the first seismic 

image of an active region on the far hemisphere of the 
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Fig. 1. — Postel projection of the far-side map of </> for 2003 
November 12 showing the strong seismic signatures close to the 
center (antipode of the Earth-facing hemisphere) of active regions 
NOAA10488(10507) and NOAA10486(10508). 

Sun (jLindsev fc Braun Il2000al) . 

The phase-sensitive holography technique is based 
on the fact that there is a phase shift (travel delay) 
between waves ente ring and exiting an active region 
(|Braun et al. I [199^ . In the analysis, the waves going 
out from a particular point in the Sun (focus) are char- 
acterized by the acoustic egression, and the waves 
going in are characterized by the acoustic ingression, 

H±{r, z,t) = j dt' j (fr' G±(r, r', z, t, t') ijj{v' , t'), 

(1) 

where r is the horizontal location of the focus, z is 
its depth, G± are Green's functions that represent the 
disturbance at (r, z, t) resulting from a unit acoustic 
impulse originating at the focus. For a detailed e xplana- 
tion of the method see lLindsev fc Braun I (|2000bD . In the 
particular case of maps of the far side, r is located at the 
surface of the non-visible hemisphere and the signal is 
measured at the Earth-facing hemisphere in an annular 
pupil centered with respect to the far-side focus. 

The far-side maps calculated from seismic holography 
represent the phase of the complex correlation between 
waves going in and out of a particular point on the non- 
visible hemisphere. In the quiet Sun, the propagation 
of the waves is well represented by the model Green's 
functions. However, when the focus is located in an 
area of high magnetic activity, the measured signal devi- 
ates from the Green's functions, effectively introducing a 
phase shift between the waves going in and out. In the 
Fourier domain, the correlation becames the product of 
the Fourier transforms, H+{r, 0, oj) and Hl{r, 0, uj) 

Civ) = I ' dio H+ir, 0, Lo)Hlir, 0, to), (2) 

and the phase (^(r) = argC(r)) is related to the per- 
turbed travel time by 

Acj) = -ljAt. (3) 

The sensitivity of the method is not uniform across the 
disk. This is due to the changing geometry of the front 
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Fig. 2. — Mean phase variation calculated from GONG far-side 
maps corresponding to the two-years interval 2003-2004. The right 
hand axis indicates the corresponding extra distance that a wave 
with a temporal frequency of 3.5 mHz will travel for the phase to 
be shifted accordingly to the left hand axis. 

side pupils as the focus move across the far-side hemi- 
sphere. Long-term averages of far-side phase maps are 
removed from individual maps to correct for this trend. 
Each far-side map is computed from a one-day series of 
1440 Global Oscillation Network (GONG) or Michelson 
Doppler Imager (MDI) Dopplergrams taken at a cadence 
of 1 minute. The Dopplergrams are taken in the photo- 
spheric line Nil A6768A. Each Dopplergram is Postel 
projected, tracking solar rotation, into a 200 x 200-pixel 
map. Waves following a ray path that bounces once at 
the solar surface before arriving to the front side (2-skip) 
are used to map the central part of the far side. A combi- 
nation of a 1-skip and 3-skip ray paths are used to extend 
the map over the full far-side hemisphere. Only waves 
with temporal frequencies between 2.5 and 4.5 mHz are 
used for the analysis. The maps are stacked into a 1440- 
frame data cube to which the analysis described above 
is applied. The resulting far-side image is itself a Postel 
projection. To reduce the errors in the calibration, only 
maps calculated from Doppler time series with a clear- 
weather duty cycle greater than 85% have been consid- 
ered. 

Figure [T] shows an example of a far-side map for 2003 
November 12 after correction. Dark areas in the map 
signify the relatively negative phase signatures, (f>, man- 
ifested by active regions. In this particular map, active 
regions NOAA 10488 and NOAA 10486 are seen on the 
far side approximately 7-8 days after they rotated across 
the West limb of the Sun, disappearing from direct view. 
They appeared on the near side in the succeeding Car- 
rington rotation. The Postel-projected far-side maps are 
then reprojected onto a longitude-sin(latitude) grid. 

Traditionally, the temporal average used to flatten the 
far-side maps is calculated for a period of low activ- 
ity. A thorough investigation of the mean phase of the 
far-side maps over long time periods has produced the 
unexpected result of a non-constant mean phase. Ba- 
sically, we find a variation of the quiet Sun from the 
model Green's functions. Figure [2] show the mean phase 
of individual farside maps calculated using Global Oscil- 
lations Network Group (GONG) data for January 2003 
to December 2004 after the average for the whole period 
has been subtracted. There is a relative variation of the 
phase of RiO.2 radians over the considered time span. A 
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Fig. 3. — Mean phase versus time. Left panel: mean phase calculated over the entire far-side hemisphere. Right panel: mean phase 
associated to high latitudes only. The right hand axis show the corresponding time delay associated with the phase shift. The change in 
the mean phase is highly anticorrelated with cycle 23 of solar activity. 



change in the mean phase can be explained by a variation 
of the cavity in which the waves propagate. 

The right-hand axis of Figure [2] shows the estimated 
distance that an acoustic wave with ^' = 3.5mHz (the 
central frequency of the wave packets considered) would 
need to travel at the photospheric sound speed (« 
8Km/s) to account for a particular phase difference. 
However, as described in the introduction, the change 
in the seismic radius is not necessarily related to a phys- 
ical change of the solar radius. It is important to note 
that the wave packets considered for calculating far-side 
maps bounce either none, one or two times before ar- 
riving on the front side and cross the subsurface layers 
several times. Hence, the expansion/contraction of the 
subsurface layers, or any other effect that produce the 
change in the seismic radius, could be amplified before 
the phase displacement is measured. 

3. SOLAR CYCLE VARIATION OF THE MEAN PHASE 

In order to search for long term changes of the mean 
phase, we examine far side maps calculated from MDI 
Dopplergrams twice daily from May 1996 to July 2008. 
An average of the complex correlation of individual maps 
calculated for all of 2007, to minimize contribution from 
activity, is used to correct for sensitivity change across 
the disk. Due to this correction, the deviation of the 
mean phase from zero is relative to the averaged period 
used for the correction, rather than from the period of 
time for which the Green's functions for the quiet Sun 
are optimal. 

Once the individual maps have been corrected, the 
mean phase is calculated for each map. Figure [3] presents 
a 30-day average of the mean phase after correction for 
the disk trend. The right axis show the corresponding 
travel time delay as per equation[3] The left panel shows 
the mean phase calculated over the full far-side hemi- 
sphere. The right panel shows the mean phase calculated 
for wave packets originating at high latitudes (> 60°). 

It is important to notice that the presence of active 
regions at the far side surface will, by itself, introduce a 
variation of the mean phase that will be correlated with 
the solar cycle. This contribution of the surface activity 
over the whole Sun is probably responsible for the differ- 



ences between the mean phase calculated using the full 
far-side hemisphere (Figure [3l left) and using exclusively 
the high latitudes (Figure[31 right), where no strong mag- 
netic activity is expected. By choosing the focus of the 
analysis to be at high latitudes, we are in effect eliminat- 
ing the in focus contribution of active regions. However, 
since the wave packets considered for far-side mapping 
bounce from one to three times at the surface on their 
way to the front side, high latitude maps are still affected 
by out of focus contribution of active regions. We ex- 
pect this contribution to be more diffuse. To accurately 
account for quiet areas only, a through analysis of he- 
lioseismic holography maps of higher resolution, where 
both the focus and the pupils are located on the front 
side, needs to be carried out. 

In order to qualitatively check the correlation of the 
mean-phase variation with the solar cycle, we define a 
Synoptic Magnetic Index (SMI) as the averaged line-of- 
sight magnetic-field strengh for individual MDI synoptic 
magnetograms. The left panel of Figure 5] shows the lon- 
gitudinally averaged line-of-sight magnetic-field strength 
for the Carrington rotations during the analyzed period 
(1909-2071). This representation shows the position and 
strength of the surface activity. We further average the 
data in latitude for each Carrington rotation to obtain 
the described SMI. Carrington rotations 1938 to 1941 
and 1944 to 1946 were eliminated of the study because 
of incomplete synoptic magnetograms. 

The right panel of Figure |4] shows a scatter plot be- 
tween the SMI and the mean phase obtained for the en- 
tire far-side hemisphere that have been averaged over 
individual Carrington rotations. The strong correlation 
is clearly seen in the figure. 

4. DISCUSSION AND CONCLUSIONS 

We find a clear anticorrelation between the seismic 
radius of the Sun and the solar activity cycle using 
a local helioseismology technique based in the analy- 
sis of wave packets. The results agre es qualita,tively 
wit h previous work by Lcfcbvrc & Koso vichev I (j2005l ) 
and K holikov fc Hill I ([20081) obtained by analysing global 
modes. 

The term seismic or acoustic radius refers to the wider 
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Fig. 4. — Left panel: Butterfly-like diagram showing the position and averaged strengh of the surface activity over time from Carrington 
rotation 1909 (May 1996) to 2071 (July 2008). Right panel: Scatter plot between the mean phase and the Synoptic Magnetic Index. The 
mean phase have been averaged over Carrington rotations for consistency with the SMI. 

the Sun and are affected by local perturbations. Hence 
the need of using local helioseismology to separate the 
relative contributions. More work is needed to disentan- 
gle the contribution of the surface magnetic areas from 
the apparent global change. Our preliminary attempt to 
isolate wave packets that have been less affected by sur- 
face activity seems to indicate an attenuation of the solar 
cycle variation in the quiet areas. Higher-resolution lo- 
cal helioseismology has the ability to separate quiet and 
active areas, which will help to further understand the 
nature of the variation. 

As a side effect of the seismic radius variation, the far 
side maps calculated from seismic holography need to 
take this temporal change into account and apply the 
appropriate correctio throughout the solar cycle in order 
to stabilize the seismic signal from active regions. 



concept of a cavity within which the waves propagate. In 
addition to a change of the physical radius, a variation 
of sound speed, temperature or changes in the superadi- 
abatic superficial layers (surface term), either global or 
localized to particular depths or latitudes, would influ- 
ence the seismic radius. Even assuming that the mean 
phase variation that we find is related to a change of the 
physical radius, th e range of change is well inside the er- 
ror bars found by iKuhn et al. I ()2004[ ) , confirming that 
the seismic inferences of the radius are more precise th an 
that of photospheric observations (jSchou et al. II1997I ). 

Because the inferred variations of the seismic solar ra- 
dius from helioseismic holography are associated with 
the propagation of wave packets, rather than individual 
modes, and the packets bounce from none to two times 
before reaching the front side where the signal is mea- 
sured, we need to be careful interpreting the results in 
terms of absolute quantities. Wave packets that prop- 
agate directly from the far-side focus to the front-side 
pupil will cross the outer layers two times, while those 
packets that bounce two times will cross the outer layers 
six times. This means our results in terms of seismic ra- 
dius variation need to be divided, on average, by a factor 
of eight, since we are comparing waves going into and out 
of the focus, resulting in inferen ces not far from those of 
global helioseismology studi es (|Lefebvre fc Kosovichev I 
l2005t iKholikov fc Hill |[200l . 

The advantage of using local-seismology methods is the 
ability to infer local properties of the Sun. The present 
work opens the road to investigate whether the varia- 
tion of the seismic radius is a global phenomenon or is 
due to variations associated with the surface magnetic 
activity. The marked difference between the mean phase 
variation from waves coming from the poles with those 
averaged over the full far-side hemisphere (see Figure [3]) 
seems to indicate that waves travelling through areas of 
high activity contribute more to the mean phase displace- 
ment, as expected. Global modes propagate throughout 
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